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We report on a lattice investigation of heavy quark diffusion within pure SU(3) plasma above the 
deconfinement transition, with the quarks treated to leading order in the heavy mass expansion. 
Using a multilevel algorithm, several volumes and lattice spacings, as well as tree-level improve- 
ment and perturbative renormalization, we measure the relevant "colour-electric" Euclidean cor- 
relator, finding that it clearly exceeds its perturbative counterpart. Even without analytic continu- 
ation, this suggests that at temperatures just above the critical one, non-perturbative interactions 
felt by the heavy quarks are stronger than within the weak-coupling expansion. After introducing 
rough modelling of the spectral shape, diffusion coefficients down to D ~ 0.5 /T appear possible. 
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Table 1: Overview of the parameter values simulated and statistics accumulated (j3 = 2N c /g 2 )- 



1. Introduction 



As was previously observed at the RHIC and recently confirmed by the heavy ion program 
at the LHC [jl|], jets containing c or b quarks (D or B mesons) get effectively "quenched", mean- 
ing that they experience a rapid kinetic equilibration with the thermal medium generated in the 
collision. The rate at which this happens appears to be much faster than at leading order of the 
weak-coupling expansion [0]. A next-to-leading order analysis suggests that indeed there are large 
corrections from higher orders [|3|], underlining the importance of a non-perturbative study. For the 
charm quark a direct measurement of the heavy quark diffusion coefficient from the current-current 
correlator is under way [Q] (with techniques developed for light quarks in ref. but given the 
systematic uncertainties involved as well as the fact that the bottom quark case is also of interest, it 
appears worthwhile to contrast this with an alternative approach, based on Heavy Quark Effective 
Theory (HQET) and valid in the large-mass limit Following a previous investigation, which 
demonstrated the principal applicability of the method [0], we report here on progress towards a 
removal of lattice artifacts and a physical interpretation of the results. 



2. Observable 



Heavy quarks carry a colour charge and, whenever there are gauge fields present, are therefore 
subject to a coloured Lorentz force. Like with other transport coefficients the corresponding "low- 
energy constants" are easiest to define at vanishing three-momentum; then the Lorentz-force is 
proportional to the electric field strength. This leads to a "colour-electric correlator" [f| 01, 



, 3 (ReTr 



U(^,x)gE i (x,0)U(z;0)gE i (0,0) 



G e (t) = --£ 

3 ;=i (ReTr[C/(i;0)] 



(2.1) 



where gE,- denotes the colour-electric field, T the temperature, and Ufa; Ti ) a Wilson line in Eu- 
clidean time direction. If the corresponding spectral function, pe, can be extracted [g], then the 
"momentum diffusion coefficient", often denoted by K, can be obtained from 

g=lim 2r PE (c a ) 

to^O CO 

According to non-relativistic linear response relations (valid for M 3> %J ', where M stands for a 
heavy quark pole mass) the corresponding "diffusion coefficient" is given by D = 2T j K. 
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Figure 1: Left: Part of a representative Monte Carlo history. Right: The corresponding averages. 

3. How to get a signal 

In order to carry out a measurement (our parameter values, referring to the standard Wilson 



gauge action, are listed in table |lj), eq. (gjj) needs to be discretized; we follow the proposal of 
ref. [§], viz- 



(3.1) 



We make use of two special techniques: the "thick" links in between the electric fields are handled 



through the link integration method [10], Ml], whereas the time intervals of width 3a, enclosing the 
electric fields, are subjected to ~ 10 extra updates with fixed boundary conditions, according to 



the multilevel philosophy [12] (a previous application at finite temperatures can be found e.g. in 
ref. Jl3|]). It depends on the parameters which of the two techniques helps more, but when both are 
combined, we always get a signal; this is illustrated in fig. [I[ In the right panel, the results have 
been normalized through [^] 
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sin 4 (7TTr) 3sin 2 (7TTr). 
4. Calibration and volume dependence 

In order to crosscheck the code, as well as for later reference, we have computed the correlator 



of eq. (3.1) to leading order in lattice perturbation theory. The result reads 

g 2 C F *[ d 3 q e *m-iT) +e qN^T ~1 
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Figure 2: Left: A crosscheck in the weak-coupling limit. Right: Volume dependence at two temperatures. 



which for N r — > °°, a — > goes over into the continuum result G^ nt (zT) of eq. ( |3.2| ). In fig. ||(left) 
a comparison of eq. ( |Q| ) with lattice measurements at very large values of j8 = 2N c /g 2 is shown, 
and we observe agreement in the weak-coupling limit. 

Proceeding towards physical measurements, fig. fright) shows results at various spatial vol- 
umes. At the current level of resolution, finite-volume effects are seen to be below statistical errors. 



5. Discretization effects 



An important systematic error originates from a finite lattice spacing (a ^ 0). A memory 
limitation currently prohibits us from increasing the spatial volume beyond 64 3 , but in view of 
the very small volume-dependence seen in fig. fright), a comparison of the lattices 48 3 x 16 and 
64 3 x 24 (the latter approximating the desired 72 3 x 24) at the same temperature allows us to test for 
the existence of scaling violations. Results at two different temperatures are shown in fig. ||(left), 
and it is clear that effects related to a finite lattice spacing need to be brought under control. 



Fortunately, the situation can be all but rectified through "tree-level improvement" [14, M 



Using eqs. (p.2|), ( |4.1| ) to determine tT from 



r LO, 
°lat 1 



(5.1) 



the pairs (zT,xT) allow us to define the tree-level improved data from the measured correlator 
Gm(zT) according to 



G imp (TT)EiG lat (TT) 



(5.2) 



The results are shown in fig. fright) and look much nicer than those in fig. ||(left). 

There are discretization effects at loop levels as well. Following the general Symanzik philoso- 
phy [16] as well as the derivation of the colour-electric correlator in ref. []6p, we may expect that the 
renormalization factor is related to the coefficient of the kinetic energy operator in lattice HQET, 
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Figure 3: Lattice spacing dependence at fixed T without (left) and with (right) tree-level improvement. 



let us denote it by c%, which needs to be appropriately tuned to match to continuum. In perturbation 
theory the value can be determined by computing the heavy quark self-energy and choosing C2 so 
as to cancel 1-loop effects specific to lattice regularization. Making use of techniques introduced 
in ref. [fnl], we find 



g 2 C F 



d 3 q 1 2 f d 4 q 1 



(27T) 3 q 2 3 J (2?r) 4 Q? 

n —n 



0.59777 



with q 2 ,Q 2 defined in eq. (4.2). The renormalization factor is then c\ « 0.83. ..0.85, i.e. quite 
modest. The results obtained after this correction (= Gi mp z) are shown in fig. @(left). Of course, 



systematic (non-perturbative) renormalization would be highly desirable [18] 



6. Physical interpretation 



The correlator Ge was computed at next-to-leading order (NLO) in continuum in ref. [19]. A 
comparison with our lattice results, improved and renormalized as explained above, is shown in 
fig. [|(left). We observe a clear enhancement over the NLO prediction over a wide Euclidean time 
interval. Subtracting the two results, one could in principle attempt a model-independent analytic 
continuation of the non-perturbative surplus, for instance along the lines of ref. [p0[]. Unfortunately, 
at the current stage, our resolution is not sufficient (~ 10 ) for this task. Therefore we resort to a 
rough model-dependent interpretation of the data in the following. 

Most transport coefficients are related to quantities which are conserved in the free limit, and 
therefore arise in connection with a narrow transport peak. This yields a contribution to the Eu- 
clidean correlator which is almost constant in T. Since the full Euclidean correlator diverges at 
short distances and is essentially perturbative there, ratios such as Gi at /GNLO should only be en- 
hanced around the middle of the Euclidean time interval (cf. e.g. refs. [|], |5]]). This is clearly not 
the case in fig. ||(left), and a narrow transport peak can be excluded, as expected ^ [T9| ]. 
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Figure 4: Left: Comparison of lattice with the NLO weak-coupling expansion JT^]. Right: Like on the left 
but normalized to the model of eq. (5.1). Closed symbols correspond to K = 4r 3 , open ones to T = 2.25T C . 



For a potentially more realistic model, we take inspiration from non-perturbative computations 
of the colour-electric correlator in other theories, one of them strongly coupled jV = 4 Super- Yang- 
Mills in the large-./V c limit 



21[ ], another classical lattice gauge theory with N c = 3 [|22|]. Both 
suggest Pe(co)/(0 to essentially flatten off below some frequency scale. Therefore, we define 

f}, 

with the free parameter K representing directly the momentum diffusion coefficient according to 



Pmodel(<f ) = max jp NL o(ft>) 



eq. ([2.2D, and compute the corresponding Euclidean correlator from 

-rr)f 



, _ rda cosh 

<Jmodel\TJ = / Pmodel(.ft>J . , m 

Jo n sinh ^ 



(6.2) 



Results are shown in fig. fright). We observe that despite extending its influence to larger fre- 
quencies than a narrow transport peak, the model does not capture the full shape of the Euclidean 
correlator; probably, extra "power" needs to be added to p((o) at even larger frequencies. 

Nevertheless, aiming at a match around the middle of the Euclidean time interval, we can 
read from fig. @(right) that K ~ 4T 3 at T = l.5T c ; K ~ 3.5T 3 at T = 2.25r c ; and K ~ 2.5T 3 at 
T = 3.0r c . The NLO weak-coupling expansion yields values k ~ 2T 3 ^ in this temperature range, 
and classical lattice gauge theory suggests that this could be an underestimate 055]]. Converting K 
to the "usual" diffusion coefficient D = 2T 2 /k, we obtain D ~ (0.5...0.8)/7\ which might lie in 
a phenomenologically acceptable range (cf. e.g. refs. [||, 23]). In ref. values D ~ 0.3/T were 
cited for the charm case. It is interesting that, despite the very rough nature of all of these estimates, 
a somewhat consistent picture appears to emerge. 



7. Conclusions 

As the results in fig. |4] show, it is possible to obtain results for the colour-electric correla- 
tor which unambiguously demand an enhancement of the non-perturbative interactions over the 
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weak-coupling description and which furthermore strongly constrain the shape of the correspond- 
ing spectral function. Of course, for quantitative statements, careful infinite- volume and continuum 
extrapolations need to be taken, statistical uncertainties should be further reduced, and a more care- 
ful modelling of the spectral shape is called for. 
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